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STUDY  OF  NONLINEAR  OSCILLATIONS  OF  ELASTIC  MEMBRANES 

1,  Foreword 

Membrane  structures  have  been  used  since  the  earliest  of  times  as  a  means  of  receiving 
acoustic  signals.  However,  in  the  past  their  analysis  relied  mostly  on  trial  and  error  while 
modem  analysis  provides  powerful  tools  for  creating  mathematical  theory  of  description 
of  their  motions.  Although  deformations  of  these  stmctures  are  essentially  nonlinear,  the 
linear  theory  provides  an  important  starting  point  for  understanding  this  complicated 
behavior. 

The  necessity  of  creation  of  passive  sensor  technologies  determines  the  importance  of  the 
study  of  elastic  membrane  oscillations.  Unattended  ground  sensors  can  have  as  their  basic 
element  a  circular  membrane  with  fixed  ends  (clamped,  simply  supported  or  damped). 
Geometry  of  the  membrane  being  chosen  to  be  circular,  the  current  investigation 
concentrates  on  the  response  analysis  of  the  membrane. 

2.  Statement  of  the  Problem,  Goals  and  the  Main  Method 

A  thin  circular  elastic  membrane  with  a  simply  supported  boundary  serves  as  a  model  of 
an  acoustic  receiver.  An  external  source  sends  a  flux  of  acoustic  waves  which  fall  on  the 
membrane  and  incite  it.  This  results  in  the  propagation  of  elastic  waves  on  the  surface  of 
the  membrane.  Vertical  oscillations  are  of  interest  and  horizontal  deflections  are 
neglected.  Nonlinear  oscillations  of  the  surface  of  the  membrane  are  governed  by  the 
damped  Boussinesq  equation.  It  is  a  second  order  in  time  and  fourth  order  in  space, 
dispersive,  dissipative,  nonlinear  equation  with  a  quadratic  nonlinearity.  Its  linear  part 
includes  the  fourth  order  classical  elasticity  operator  and  the  third  order  term  responsible 
for  internal  friction.  Quadratic  nonlinearity  present  in  the  equation  includes  second 
derivatives.  It  reflects  dependence  of  the  vertical  deflection  on  the  local  curvature  of  the 
membrane. 

Goals  of  investigation  consist  in  developing  an  efficient  algorithm  of  solving  the  direct 
nonlinear  boundary-value  problem  in  question,  i.e.  deriving  the  formulas  for  the  vertical 
deflection  of  the  membrane,  proving  the  convergence  of  the  corresponding  series  in  a 
certain  function  space  and  finding  the  acoustic  pressure.  Conducting  numerical 
simulations  is  necessary  for  the  comparison  of  behavior  of  the  linear  and  nonlinear 
membranes.  In  addition,  a  progress  should  be  made  in  finding  the  direction  to  the 
acoustic  source  and  exploiting  the  nonlinearity  for  making  the  membrane  more  sensitive. 
The  latter  goal  stands  out  as  a  new  type  of  an  inverse  nonlinear  problem  which  has  never 
been  solved  before. 

Nonlinear  version  of  the  method  of  eigenfunction  expansions  is  chosen  as  a  main  tool  of 
constructing  solutions.  This  method  is  well  known  for  the  classical  linear  models  of  wave 
propagation.  Extending  the  method  for  solving  nonlinear  boundary  value  problems  has 


been  a  reeent  development  made  by  the  Prineipal  Investigator  of  the  present  work. 
Solution  is  obtained  in  the  form  of  the  series  of  eigenfunetions  of  the  Laplaee  operator  in 
a  dise  with  homogeneous  boundary  eonditions  and  periodieity  eonditions  in  the  angle. 
More  preeisely,  it  is  a  Fourier-Bessel  series,  where  Bessel  funetions  are  the  radial 
eigenfunetions  and  exponentials  are  the  angular  ones.  Bessel  funetion  zeros  depend  on 
two  indiees  whieh  refleets  the  interaetion  of  radial  and  angular  eigenfunetions.  Expanding 
the  nonlinearity  into  the  eigenfunetion  series  and  obtaining  the  estimates  of  the 
eigenfunetion  eoeffieients  is  the  key  issue  in  applying  this  method.  In  eontrast  to 
Galerkin’s  method  whieh  allows  only  to  prove  existenee,  our  approaeh  allows  to 
eonstruet  solutions. 


3,  Results 

Proving  eonvergenee  of  the  series  representing  the  solution  and  the  nonlinearity  forms  the 
main  diffieulty  and  defines  the  funetion  spaee  where  the  solution  exists.  It  was 
established  that  elassieal  solutions  of  the  boundary  value  problem  ean  not  be  eonstrueted 
in  prineiple  sinee  eonvergenee  of  the  series  is  not  enough  to  provide  the  neeessary 
regularity.  This  is  typieal  for  many  nonlinear  initial-boundary-value  problems.  However, 
mild  solutions  ean  be  eonstrueted.  They  are  defined  as  solutions  of  the  integral  equation 
obtained  by  integrating  the  original  boundary-value  problem  with  respeet  to  time. 

New  funetion  spaees  were  introdueed,  namely  anisotropie  Sobolev  spaees.  The  norm  in 
sueh  a  spaee  is  the  sum  of  the  L2-norms  of  the  usual  derivatives  weighted  by  tangential 
derivatives.  These  spaees  are  sensitive  to  the  energy  transfer  in  the  angular  direetion  and 
refieet  the  interaetion  of  angular  waves  and  the  radial  ones. 

A  new  family  of  speeial  funetions  was  introdueed  for  improving  the  eonvergenee  of  the 
series  representing  the  solution.  They  owe  their  appearanee  to  the  nonlinearity  and  the 
eireular  geometry.  These  funetions  appear  as  a  result  of  the  series  multiplieation.  They 
reeeived  the  name  of  eonvolutions  of  Rayleigh  funetions  with  respeet  to  the  Bessel 
index.  Rayleigh  funetions  are  well  known  in  the  elassieal  theory  of  linear  dise 
oscillations.  They  are  defined  as  series  of  inverse  powers  of  the  Bessel  function  zeros  and 
reflect  the  circular  geometry  for  linear  problems  and  the  boundary  conditions. 
Convolutions  of  these  functions  reflect  circular  geometry,  boundary  conditions  and  the 
nonlinearity.  These  functions  allow  to  improve  the  convergence  of  the  eigenfunction 
series  representing  the  solutions  and  to  improve  greatly  the  computational  effectiveness 
of  the  algorithm.  Convolutions  of  Rayleigh  functions  were  expressed  in  terms  of  the  psi- 
function  (logarithmic  derivative  of  the  gamma- function  and  its  derivatives).  Asymptotic 
expansions  for  large  values  of  the  angular  index  were  obtained.  These  asymptotics  are 
essential  for  the  estimates  of  the  Fourier-Bessel  coefficients  in  the  representations  of 
solutions. 

Eigenfunciton  series  representations  were  obtained  for  the  vertical  deflection  of  the 
circular  membrane  and  the  acoustic  field  in  the  plane  of  the  membrane.  These  are  mild 
solutions  of  the  nonlinear  boundary  value  problem  in  question.  If  a  typical  time-harmonic 


signal  is  chosen  as  a  source,  then  the  final  formulas  show  a  typically  nonlinear  effect  of 
frequency  multiplication. 

Numerical  simulations  were  conducted  in  order  to  compare  linear  and  nonlinear 
oscillations.  Truncated  eigenfunction  series  solutions  were  used  for  computations  in  the 
nonlinear  case.  Different  angles  for  incident  acoustic  waves  were  tried.  It  was  established 
that  the  created  algorithm  has  excellent  convergence  properties.  Just  a  few  iterations  were 
needed  for  securing  enough  precision.  The  conducted  simulations  clearly  showed  the 
picture  of  nonlinear  waves  traveling  on  the  surface  of  the  membrane  and  a  much  better 
sensitivity  of  the  nonlinear  membrane  in  comparison  with  the  linear  one. 

As  regards  the  most  difficult  problem  of  finding  the  angle  to  the  acoustic  source  (inverse 
nonlinear  problem  which  deserves  the  name  of  the  problem  of  the  century),  some 
progress  has  been  made.  Indeed,  one  can  establish  the  plane  where  the  source  lies.  On  the 
border  (perimeter)  of  the  membrane  there  is  a  spot  where  the  maximum  of  the  wave 
energy  is  located.  One  can  draw  the  diameter  line  through  the  center  of  the  disc  and  the 
point  of  the  maximum.  The  plane  going  through  this  line  perpendicular  to  the  membrane 
contains  the  source  of  acoustic  waves.  Rigorous  mathematical  proof  of  this  fact  is  being 
developed. 
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